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ABSTRACT 

We consider particle acceleration by large-scale incompressible turbulence with a lengthscale larger 
than the particle mean free path. We derive an ensemble-averaged transport equation of energetic 
charged particles from an extended transport equation which contains the shear acceleration. The 
ensemble-averaged transport equation describes particle acceleration by incompressible turbulence 
(turbulent shear acceleration) . We find that for Kolmogorov turbulence, the turbulent shear acceler- 
ation becomes important in small scale. Moreover, by Monte Carlo simulations, we confirm that the 
ensemble-averaged transport equation describes the turbulent shear acceleration. 
Subject headings: acceleration of particles — turbulence — plasmas — cosmic rays — ISM: supernova 
remnants 
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1. INTRODUCTION 

Charged particles are accelerated to relativistic 
energies in many astrophysical objects. In ad- 
dition, turbulence is also expected. In fact, 
strong turbulence is observed in recent two and 
three di mensional simulation s for supernova remnants 
/SNRs) (IGiacalone fc Jokbii '200?'; 'In oue et all [200l 
ICuo et all 120121: ICaorioh & Spitkovskv Tl2nTll. pul- 
sar wi nd nebulae (PWNe) pC pmissarov & Lvubarskv ' 
120041 : iDel Zanna et al.l l2004rTPorth et al. 2012.\ . as- 



trophysicai iets (lAlovet all 119991: iMizuta et "ail 120101: 
iLopez-Camara et al.ll2012[ ). etc. 

There are mainly two acceleration mechanisms by 
turbulence. One is due to wave-particle interac- 
tions, where the particle mean free path is compara- 
ble to the wavelength of electromagnet i c fluc tuations 
(e.g ISkiUind [TQTI ISchlickeiser fc Milled [l998l) . The 
other is due to large-scale fluctuations of plasma flows, 
where the particle mean free path is small er than 
turbulent scales (e.g. Bv kov fc Toptvginl 119931 ). Tur- 
bulence is generally divided by compressible and in- 
compressible modes. Particle acceleration by large- 
scale compressib le turbulence has been discussed in - 
many paper (e.g. jBvkov fc Toptvgin|[l982tlPtuskinlll988l: 
iJokipii fc Led [2010f r However, particle acceleration by 
large-scale incompressible turbulence (turbulent shear 
acceleration) has not bee n investigated in detail, while 
iBvkov fc Toptyginl (|1983[ ) has briefly discussed the tur- 
bulent shear acceleration. 

Particle acceleration by a simple incompressible flow 
(shear flow) has already investigated in many paper (e.g. 
Berezhko fc Krvmskii 1981; Earl et al. 1988; Webb 1989; 
Ostrowskil [19901 : iHieger fc Dufrvll200Q) . However, shear 



flows are potentially unstable to the Kelvin-Helmholtz 
instability and produce turbulence. Therefore, the tur- 
bulent shear acceleration is expected to be important. 
In this Letter, we investigate the turbulent shear accel- 
eration by considering ensemble average of an extended 
transport equation which includes particle acceleration 
by shear flows. 

We first derive an ensemble-averaged transport equa- 
tion in Section 2, and provide its analytical solutions for 



simple cases in Section 3. We then perform Monte Carlo 
simulations in Section 4. Section 5 is devoted to the dis- 



2. DERIVATION OF THE ENSEMBLE-AVERAGED 
TRANSPORT EQUATION 

In this section, we derive the ensemble-averaged trans- 
port equation of energetic particles. Propagation and ac- 
celeration of energetic charged particl es in a p lasma flow 
are described by a transport equation. iParker I (1965) de- 
rived the transport equation which includes spatial diffu- 
sion, convection, and adiabatic acceleration. After that 
his work was extended by several authors. For isotropic 
diffusion and a nonrelativistic plasma flow, the extended 
transp ort equation is given by (Equation (4.5) o f lWebbI 
(fT989h and Equation (9) of lWilhams eTall (fl99l ) 
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F{p, X, t), Xi, Ui, k(p) are the distribution function, posi- 
tion, plasma velocity, and spatial diffusion coefficient, re- 
spectively. V and p are the particle velocity and four mo- 
mentum in the fiuid rest frame, respectively. The spatial 
diffusion coefficient, k, is represented by k(j>) = t(p)w^/3 
for isotropic diffusion, where t{p) is the mean scattering 
time and tv is the particle mean free path. The first 
four terms of Equation ([ij are the same as the Parker 
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equation and the others are additional terms. The fifth 
term describes the shear acceleration and the sixth term 
becomes important for v ^ Ui. 

In order to understand essential features of the turbu- 
lent shear acceleration, we consider incompressible tur- 
bulence (dUi/dxi = 0) and do not take into account the 
spatial transport, that is, we consider the spatially aver- 
aged distribution function, Jy Fd^x, where V is the 
system volume that we consider. By integrating Equa- 
tion ([T]), the extended transport equation can be reduced 
to 
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where Q{p) is the particle flux passing through the sur- 
face of the integrated volume. As long as we consider 
a timescale smaller than V^^'^/Ui and V^^^/k, we can 
neglect the particle flux, Q{p). In other words, we can 
neglect escape of particles from the system when we con- 
sider a sufficiently large system size. 

In this Letter, we assume that the plasma velocity field, 
Ui{x,t), is static, random, statistically homogenous and 
isotropic incompressible turbulence, that is, Ui = 5ui(x) 
and {6ui) — 0, where (...) denotes ensemble average. The 
correlation function of the plasma velocity field is given 
by 

(5«,(x)5^,(x')) = I ^X.,(k)e^{'='(---;) , (4) 



and 



(5) 



where k and S{k) are the wavenumber and spectrum of 
incompressible turbulence, respectively. As long as we 
consider only particle acceleration, we can assume the 
velocity field to be static when the scattering timescale, 
r, is smaller than the variable timescale of fluid, T ^ 
{k X max Uph})~^, where Uph is a phase velocity. In 
this letter, we consider rvk < 1 and v > max Wph}, 
so that t/T ^ rk X max {Su, Wph} < Tvk < 1. Hence, we 
can assume a static velocity field in this letter. 

The distribution function of particles can also be di- 
vided by an ensemble-averaged component and a fiuctu- 
ated one, that is, F = f + Sf and (F) = f. The spatial 
average in Equation ([3]) can be interpreted as ensemble 
average because we consider a system size larger than the 
turbulent scale. Then, from Equation ([3]), the ensemble- 
averaged transport equation is represented by 
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where we have assumed that distributions of 6f and 6ui 
are symmetric about the mean values, / and 0, respec- 
tively, so that third moments are zero. From Equations 
([4]), ([5]), and dH), the ensemble- averaged transport equa- 



tion can be represented by 
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where the momentum diffusion coefficient, Z^tsaIp), is 
given by 

D.sAip) = -^pMp) J J^Sik)k^ [^-+'-^ 

(8) 

We here consider turbulence with a large lengthscale 
compared with the particle mean free path, rw, so 
that the upper limit of fc-integral should be limited by 
minjfcniax, fcrcs} whcrc fcmax is thc maximum wavenum- 
ber of turbulence and fcics ~ {'tv)^^ ■ The momentum 
diffusion coefficient, D^sA, is dominated by small scale 
turbulence when k^S{k) is an increasing function of k. 
Therefore, the turbulent shear acceleration becomes im- 
portant in the small scale for a Kolmogorov-like spectrum 
{S{k) oc fc^"/3). 

3. ANALYTICAL SOLUTION 

In this section, we present specific expressions of the 
momentum diffusion coefficient and analytical solutions 
of the ensemble-averaged transport equation for simple 
velocity spectra. We especially focus on the turbulent 
shear acceleration of relativistic particles (u « c) in non- 
relativistic turbulence {{Su'^) ^ c^), so that we neglect 
the term of {Su^)/v'^ in Equation We here assume a 
functional form of the mean scattering time, t{p), to be 
tq{p/po)°' , where po and tq are the initial four momen- 
tum and the mean scattering time of particles with pq, 
respectively. To make the expression simple, hereafter 
the four momentum, time, and momentum diffusion co- 
efficient are normalized by Po,to, and p^r^^, respectively. 
Normalized quantities are denoted with a tilde. 

For a static monochromatic spectrum of incompressible 
turbulence, S{k) is given by 
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From Equations ([5]) and the momentum diffusion 
coefficient is represented by 
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For a static Kolmogorov-like spectrum of incompress- 
ible turbulence, we assume that S{k) is given by 



((5w2)(27r)2 
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fc-"/3 (for ko<k< fc,nax) 

(11) 



Then, from Equations ([5]) and dTll . the momentum dif- 
fusion coefficient is represented by 



D 



TSA 



{Tavkof {Su^) ( min{fcmax,fcros}V^^ 
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(12) 

where we have assumed fco ^ min{fci„ax, fcres}- The 
factor, (min{fcmax, fcres}/fco)''''^j is expected to be large. 
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Therefore, the Kolmogorov-like turbulent cascade en- 
hances the turbulent shear acceleration. For fcrcs < fcmax, 
-Dtsa is represented by 



D. 
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30 
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(13) 



Therefore, the momentum diffusion coefficient can be 
represented by -Dtsa — D^p^^^ for above simple cases, 
where /3 = a for the monochromatic spectrum and the 
Kolmogorov spectrum of the case A^rcs ^max; and 
/3 = —a/3 for the Kolmogorov spectrum of the case 

^rcs ^ ^max- 

We next discuss analytical solutions of the ensemble- 
averaged transport equation. We assume that particles 
are uniformly distributed in the three dimensional space 
and injected at time, t — 0, with the four momentum, 
p = 1. Then, the ensemble-averaged transport equation 
is represented by 
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S{i)S{p~l) , (14) 



where N is the number of injected particles. If the mo- 
mentum diffusion coefficient is represented by Dtsa = 
DqP^^^, for /3 ^ 0, the so lution is given by ([Berezhkol 
[1981 IRieeer fc Duffvl[2006l) 



N 



4,TT\P\Dot 

X/|l+3//3| 



^,p-(3+«/2exp 



lj-pp_ 

13^ Dot 



(15) 



where is the modified Bessel function of the first kind. 
The solution approches p^f oc jp"^ for jp > 1. For /? = 0, 
the solution is given by (|Rieger fc Duffvl[2nM ) 
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(16) 



and the evolution of the mean momentum, Pmit) 
N^^ J pf{p,t)4:TTp^dp, is given by 



Pm(t) = exp 



(17) 



Note that solutions of Equations ()15|) and ([161) ^-^^ 
valid for t ^ 1 and p ^ 1 because of causality. 

4. MONTE CARLO SIMULATION 

In order to confirm analytical solutions presented in the 
previous section, we perform test particle Monte Carlo 
simulations. We here focus on static, statistically ho- 
mogenous and isotropic incompressible turbulence, that 
is, the velocity field, (5Mi(x), is divergence free. Such a 
vector field is numerically constr ucted by a summation of 
many transverse waves ( Giaca lone fc Jokip ii 1999). Sim- 
ulation particles are isotropically and elastically scat- 
tered in the local fluid frame and move in a straight line 
between each scattering. The mean scattering time is 
given by r = To{p/po)°'- We use 10'' simulation par- 
ticles with the initial four momentum po = lOmc and 




Fig. 1. — Time evolution of the mean four momentum for a = 
and TQcko = 10~^. The dots and solid lines show the resul ts o f 
Monte Ca rlo simulations and analytical solutions ( Equations ljlO|l , 
II12I I and I I17I I). respectively. The red and blue show cases of the 

monochromatic and Kolmogorov spectra with rocfcmax = 10"^'^^, 
respectively. 




TnCkn 



Fig. 2. — Wavenumber dependence of the momentum diffusion 
coefficient for a = 0. The dots and solid lines show the results 
of M onte Car lo simulations and analytical solutions of Equations 
moi l and I I12I I. respectively. The red and blue show cases of the 
monochromatic and Kolmogorov spectra with Tocfcmax = 10"^''^, 
respectively. 



10'' 




Fig. 3. — Distribution function at t/ro = 10"" for q = and 
TQcko = 10-1. The histo grams and solid lines show the results of 
Monte Carlo simulations and analytical solutions of Equations ||16|I , 
respectively. The red and blue show cases of the monochromatic 
and Kolmogorov spectra with Tocfcmax = IQ— respectively. 
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Fig. 4. — Distribution function for the monochromatic spectrum 
with TQcko = 0.01 and a = I. The red histograms and black solid 
lines show the results of M onte Carlo simulations and analytical 
solutions of Equations II15I I. respectively. The thin and thick show 
distribution functions at t/ro = 5 X 10® and lO'', respectively. 

100 transverse waves in order to construct velocity fields, 
where m and c are the particle mass and the speed of 
light. The mean amplitude of velocity fluctuations is 
taken to be (Su^) = (0.05c)^. We set the maximum 
wavenumber to be Tocfcmax = 10"^^"^ for the Kolmogorov 
spectrum. 

We first discuss results of Monte Carlo simulations for 
the momentum-independent scattering, that is, a — /3 = 
0. Figure 1 shows the evolution of the mean four momen- 
tum for a = and Tgcfco = 10^^. Particles are acceler- 
ated, and simulation results are in good agreement with 
analytical solutions of Equations (ITUl) , (|12l) and ([T7| . By 
comparing the growth rate of the mean momentum of 
simulation particles with Equation (1171) . we can obtain 
the momentum diffusion coefficient of Monte Carlo sim- 
ulations. 

Figure 2 shows the wavenumber dependence of the mo- 
mentum diffusion coefficient, Dq = Dtsatq / for a = 0. 
Simulation results are in good agreement with Equations 
pH]) and as long as tqc/co < 1, but simulation results 
for the monochromatic spectrum deviate from Equation 
(fT0|) at Tgcfco > 1. As already mentioned in Section 2, 
this is because our treatment is not valid when the parti- 
cle mean free path is larger than the turbulent scale. Fur- 
thermore, we have confirmed that the Kolmogorov-like 
turbulent cascade (blue) enhances the turbulent shear 
acceleration. 

Figure 3 shows the distribution function, pdN/dp oc 
p^fip, t), at t/TQ = 10^ for a = and rocfco = 10^^ Sim- 
ulation results (histograms) are in excellent agreement 
with analytical solutions of Equation ([16]) (solid lines) for 
the monochromatic (red) and Kolmogorov (blue) spec- 
tra. 

Figure 4 shows the distribution function at I/tq — 
5 X 10^ and 10^ for the monochromatic spectrum with 
TQcko = 10^^ and the Bohm-like diffusion, that is, 
a — P = 1. Simulation results (histograms) are in ex- 
cellent agreement with analytical solutions of Equation 
(solid lines) except for above p/po ^ 10^. As men- 
tioned above, the disagreement is due to T{p)vkQ > 1 at 

p/po > 10^. 

Therefore, we have confirmed that the ensemble- 
averaged transport equation for incompressible turbu- 
lence describes the turbulent shear acceleration and, that 



is valid as long as the turbulent scale, ~ fc ^, is larger 
than the particle mean free path, tv. 

5. DISCUSSION 

We first discuss another importa nt effect of turbu- 
lence on the particle transport. IBvkov fc ToptvginI 
(jl993[) shows that turbulence enhances spatial diffusion. 
For strong turbulence, Kturb becomes of the order of 
Lo^/iS^ (jBvkov fc ToptveinllTOOl ). where Lq is the in- 
jection lengthscale of turbulence, so that spatial diffusion 
of particles with a small mean free path is dominated by 
turbulent diffusion and an energy-independent diffusion 
is realized. The ratio of the turbulent diffusion and the 
Bohm diffusion, Ksohm, is given by 



^^^3x10^ 



'^Bohm 



m-oC 



B 



\ c / V 1 / \ 1 pc 

(18) 

where Wp and B are the proton mass and magnetic field, 
respectively. Therefore, turbulent diffusion of energetic 
particles could be important in SNRs, PWNe, astrophys- 
ical jets, etc. 

From Equation ()13|) . the acceleration timescale, iacc = 
P^/Dtsa, of the turbulent shear acceleration for the Kol- 
mogorov spectrum of the case fcjcs < ^max is represented 

by 

,2 



30 



= 9 X 10^ sec 

p Y^^'fiSu^) 



t{p) 



B 



1 AiC 



-1/3 / ^ X 2/3 



1 pc 



(19) 



where we have assumed Lq = 27r/fco and the Bohm dif- 
fusion, t{p) = p/(eB), in the last equation. Therefore, 
particles can be accelerated to relativistic energies by 
large-scale turbulence in many astrophysical objects. In 
addition, if particles are initially accelerated at the shock, 
large-scale turbulence can change energy spectra of the 
accelerated particles in the shock downstream region. 

Next, we compare the turbulent shear acceleration and 
particles acceleration by small-scale incompressible tur- 
bulenc e, that is. the s econd order acceleration by Alfven 
waves (jSkillindll975D . The momentum diffusion coeffi- 
cient of the second order acceleration by Alfven waves is 
given by Da ^ p'^v'a/ (9^), where va is the Alfven veloc- 
ity. The ratio of the turbulent shear acceleration and the 
second order acceleration by Alfven waves is given by 



-Ptsa (Sv^ 
Da ^ vl 



2/3 



(20) 



where we have adopted Equation (fTn|) as -Dtsa- There- 
fore, the turbulent shear acceleration could be more effi- 
cient than the second order acceleration by Alfven waves 
for super- Alfvenic turbulence {^/{Sii^ > va(tv/ Lq)^/^) 
. In other words, the turbulent shear acceleration be- 
comes important when there are strong magnetic field 
fluctuations {SB/Bq > 1) because the plasma veloc- 
ity fluctuation by Alfven waves, Su, is represented by 
Su = va{SB / Bq), where 6B and Bq are the fluctuated 
and mean magnetic fields, respectively. Such a situation 
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is expected to be realized in the d ownstream region of 
a high Alfven Mach nu mber shock ([Giacalone fc Jokipiil 
[20071 llnoue erall[2009l) . 

We have considered only isotropic diffusion and non- 
relativistic turbulence in this Letter. Spatial diffu- 
sion is generally anisotropic because of the magnetic 
field. Isotropic diffusion is realized when magnetic 
field fluctuations with lengthscale comparable to the 
particle mean free path a re large {SB/Bq > I) (e.g. 
iGiacalone fc Jokipiil Il999t ). Therefore, as discussed 
above, the turbulent shear acceleration is important 
when isotropic diffusion is realized. Simple extensions 
to anisotropic diffusion and relativistic turbulence are 
straightforward because e xtensio n s of E q uation ([H) have 
already been provided by iWebbI ()1989D : iWilliams et al.l 
()1993[ ). This calculation will be addressed in the future 
work. 

6. SUMMARY 

In this Letter, we have derived a particle transport 
equation averaged over random plasma flows in order to 



understand particle acceleration in incompressible tur- 
bulence with a larger lengthscale than the particle mean 
free path. We have considered ensemble average of the 
extended trans p ort eq uation provided by Webb (198£i); 
IWilhams et all (fl993h . This is a simple extension of 
previous work that considered ensemble aver age of the 
transport equation provided by iParkeTI (|1965[ ). We have 
found that the turbulent shear acceleration by incom- 
pressible turbulence becomes important in small scale 
for Kolmogorov-like turbulence. Moreover, we have per- 
formed Monte Carlo simulations and confirmed the tur- 
bulent shear acceleration. Recent simulations show that 
turbulence is produced in many astrophysical objects, so 
that turbulent diffusion and turbulent acceleration are 
expected to be important. 
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Education, Culture, Sports, Science, and Technology 
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